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2^ ; Abstract 

In the framework of the simplest little Higgs model (SLHM), we study the production of a pair 
of neutral CP-even Higgs bosons at the LHC. First, we examine the production rate and find that 
it can be significantly larger than the SM prediction. Then we investigate the decays of the Higgs- 
pair and find that for a low Higgs mass their dominant decay mode is hh — > nnnn (n is a CP-odd 

(N ■ 

scalar) while hh — > bbnn and hh — > nnWW may also have sizable ratios. Finally, we comparatively 

i> : _ 

fNJ , study the rates of pp — > hh — > bbr + r , pp — > hh — > 6677, and pp — > hh — > WWWW in the 

OC • 

SLHM and the littlest Higgs models (LHT). We find that for a light Higgs, compared with the SM 

00 ■ 

predictions, all the three rates can be sizably enhanced in the LHT but severely suppressed in the 

On i 

SLHM; while for an intermediately heavy Higgs, both the LHT and SLHM can enhance sizably 

£> 

' the SM predictions. 
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I. INTRODUCTION 

Little Higgs theory |l| has been proposed as an interesting solution to the hierarchy prob- 
em. So far various realizations of the little Higgs symmetry structure have been proposed 
, |5j, which can be categorized generally into two classes [6|. One class use the product 
group, represented by the littlest Higgs model [4j, in which the SM SU(2)l gauge group is 
from the diagonal breaking of two (or more) gauge groups. The other class use the simple 
group, represented by the simplest little Higgs model (SLHM) Q], in which a single larger 
gauge group is broken down to the SM SU(2)l. Of course, different realizations give different 
phenomenology, which will be tested at the LHC 

Since these little Higgs models mainly alter the properties of the Higgs boson and the 
top quark hints of these models may he unravelled from various H.ggs hosou and top quark 
processes |7j. The Higg s-pair production at the LHC, albeit with a small production rate, 
is rather important because it will provide a way to probe the Higgs self-coupling A. With 
the designed luminosity, it is possible for the LHC to establish that the SM Higgs boson has 
a non-zero self-coupling and the ratio X/Xsm can be restricted to a range of — 3.7 at 95% 
confidence level if its mass is between 150 GeV and 200 GeV gj]. Such Higgs-pair production 
is sensitive to new physics and has been studied in various new physics models 9J. In the 



littlest Higgs models without and with T-parity, this process was studied in |10J and 11], 
respectively. In this work, we study this process in the SLHM. We will first examine the 
Higg-pair production rate in the SLHM and compare with the SM prediction. Then we 
study the decays of the Higgs-pair. Finally, we study the rates of pp — > hh — » bbr + T~ (6677) 
and pp —>■ hh — > WWWW , comparing the prediction of the SLHM with the littlest Higgs 
models. 

This work is organized as follows. In Sec. II we recapitulate the SLHM. In Sec. Ill we 
calculate the Higgs-pair production cross section at the LHC. In Sec. IV, we study the decays 
of the Higgs-pair and the rates of pp —>■ hh — > bbr + r~ (bbry-f) and pp — > hh —>■ WWWW . 
Finally, we give our conclusion in Sec. V. 
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II. SIMPLEST LITTLE HIGGS MODEL 



The SLHM is based on [SU (3) x U{l)x} 2 global symmetry. The gauge symmetry SU (3) x 
U(l)x is broken down to the SM electroweak gauge group by two copies of scalar fields $1 
and $ 2 , which are triplets under the SU(3) with aligned VEVs f\ and / 2 . The uneaten five 
pseudo-Goldstone bosons can be parameterized as 



where 
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/ = V fx + fi an d tp = tan/3 = f^/ fx- Under the SU(2) L SM gauge group, n is CP-odd 
singlet, while H transforms as a doublet and can be identified as the SM Higgs doublet. The 
kinetic term in the non-linear sigma model is 
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5=1,2 



(3) 



where ^ = gtanOw/y 1 — tan 2 0yj//3 with being the electroweak mixing angle. As $i 
and $2 develop their VEVs, the new heavy gauge bosons Z', Y°, and W ± get their masses 
proportional to /. 

The gauged 577(3) symmetry promotes the SM fermion doublets into SU(3) triplets. 
There are two possible gauge charge assignments for the fermions: the 'universal' embed- 
ding and the 'anomaly-free' embedding. The first choice is not favored by the electroweak 
precision data 0, so we focus on the second way of embedding. The quark Yukawa inter- 
actions for the third generation and the first two generations can be written respectively 
as 



c 3 



i\{t c MQ 3 + i\\t c 2 & 2 Qz + i'-^d c m e ijk <S>\& 2 Ql + h.c. 
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where n = 1,2 are the first two generations indices; i,j, k = 1,2,3; Q3 = {ti^L^T^} and 
Qn = {dnL, —UnL, iDni}] d c m runs over (d c , s c , b c , D c , S c ); d\ n and d\ n are linear combinations 
of d c and D c for n = 1 and of s c and S c for n = 2; runs over (m c , c c , t c , T c ). For simplicity, 
we assume the quark flavor mixing are small and neglect the mixing effects. From Eqs. (T4|) 
and ([5]), we can get the Higgs boson interactions and the mass terms for the three generations 
of quarks: 

C t ~ -/A* [xlc p tl{-s x t L + Cl T L ) + s p t c 2 (s 2 t L + c 2 T L )] + /i.e., (6) 
C dn ~ -/A* 1 [x^c^Csid,^ + cxAa) + s p d c 2 (-s 2 d nL + c 2J D nL )] + /i.e., (7) 
£ g ~ —^-f 2 Sf3Cf3S 3 q c q L + h.c. (q = u, c, 6) (8) 



where 



t _ M d n _ ^1 _ /g _ /l 



. t p {h + v) . (h + v) . (h + v)(t 2 + l) 

S] sm ^r> aassu w 53 " 8111 V2W ' (9) 

with h and t> being the neutral Higgs boson field and its VEV, respectively. The mass 
eigenstates are obtained by mixing the corresponding interaction eigenstates, e.g., the mass 
eigenstates (t mL ,T mL ) and (t^,T^) are respectively the mixtures of (tL,T L ) and (t c ,T c ). 
The diagonalization of the mass matrix in Eqs. ([6]) and ([7j) was performed numerically in our 
analysis, and the relevant couplings with Higgs boson can also be obtained without resort 
to any expansion of v/f. Hereafter we denote the mass eigenstates without the subscript 
''m 1 for simplicity. 

The Yukawa and gauge interactions break the global symmetry and then provide a poten- 
tial for the Higgs boson. However, the Coleman- Weinberg potential alone is not sufficient 
since the generated Higgs mass is too heavy and the new CP-odd scalar rj is massless. 
Therefore, one can introduce a tree-level \i term which can partially cancel the Higgs mass 

" ^ + = "V/V.COS (jg-.) COS . (10 ) 

The scalar potential becomes 

V = -m 2 H ] H + \(H ] H) 2 - ^m 2 r] 2 + X'H^Hrj 2 + ■■■ , (11) 
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where 



2 2 A* 

m =m , A = Ao - — 3 3 , 

^s%c%f 2 



A' 



(12) 



with mo and Ao being respectively the one-loop contributions to the Higgs boson mass and 
the quartic couplings from the contributions of fermion loops and gauge boson loops Q]. 
The Higgs VEV, the Higgs boson mass and the mass of r\ are given by 



m 



m h = 2m , mf 



cos 



The Coleman- Weinberg potential involves the following parameters: 



(13) 



/, x\, tp, fi, m v , m h ,v. 



Due to the modification of the observed W gauge boson mass, v is defined as 
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(15) 



where Vq = 246 GeV is the SM Higgs VEV. Assuming that there are no large direct contri- 
butions to the potential from physics at the cutoff, we can determine other parameters in 
Eq. (|14l) from /, tp and m^ with the definition of v in Eq. ffTBI) . 



III. HIGGS-PAIR PRODUCTION AT LHC 

At the LHC the Higgs-pair production can proceed through gluon-gluon fusion and bb 
annihilation, as shown in Figs. (pQ) and ([2]), respectively. For the bb annihilation process, 
the SLHM can give the additional contributions through the tree-level hhbb coupling and 
the modified hbb coupling. For the gluon-gluon fusion process, the top-quark loops give 
additional contributions through the tree-level hhtt coupling and the modified hit coupling. 
In addition to the top-quark loops, the loops of the new heavy partner quarks T, D and 
S also come into play. Due to the large top quark Yukawa coupling and the large parton 
distribution function of gluon at the LHC, the contributions of the gluon-gluon fusion process 
can be dominant over bb annihilation process. 

ig. Q] are straightforward. Each loop diagram 



The calculations of the loop diagrams in 
is composed of some scalar loop functions 

[ijj]. The calculations are tedious and the analytical expressions are lengthy, which are not 



13] which are calculated by using LoopTools 
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FIG. 1: Feynman diagrams for Higgs-pair production via gluon-gluon fusion in the SLHM. Here 
i,j = 1,2 with (/i,/2) denoting (t,T) or (d,D) or (s,S). The diagrams by exchanging the two 
gluons or exchanging the two Higgs bosons in (c,d) are not shown here. 



(a) 



(b) 





(c) (d) 
FIG. 2: Feynman diagrams for Higgs-pair production via bb annihilation in the SLHM. 

presented here. The hadronic cross section at the LHC is obtained by convoluting the parton 
cross section with the parton distribution functions. In our calculations we use CTEQ6L 2] 
to generate the parton distributions with the renormalization scale /j,r and the factorization 
scale llf chosen to be fiR = llf — 2rrih and use the two-loop running coupling constant a s 
with a 3 {mz) = 0.118. 

The SM input parameters relevant in our study are taken as m t = 171.2 GeV and 
mz = 91.1876 GeV [16]. The free SLHM parameters are f,tf3,rrih,x'l('mo) and x\(ms)- As 
shown above, the parameters x l Xl li, m v can be determined by /, tp, rrih and v. The small 
mass of the d (s) quark requires one of the couplings Xf and (AJ and A|) to be very small, 
so there is almost no mixing between the SM down-type quarks and their heavy partners. 
We assume \f (A*) is small, and take x x = 1.1 x 10 -4 and x s x = 2.1 x 1CT 3 , which can 
make the masses of D and S in the range of 1-2 TeV with other parameters fixed as in our 
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FIG. 3: Hadronic cross section of Higgs-pair production at the LHC versus the Higgs boson mass. 

calculations. In fact, our results show that the contributions from d and D (s and S) are 
small compared with the effects from t and T. So, different choices of xf and x s x do not have 
sizable effects on our results. 

Electroweak precision data can give the strong constraints on the scale /. The |5| shows 
that the LEP-II data requires / > 2 TeV. In addition, the contributions to electroweak 
precision data can be suppressed by large tp. Ref. [17| gives a lower bound of / > 4.5 
TeV from the oblique parameter S while a recent study of Z leptonic decay gives a stronger 



bound of / > 5.6 TeV [18]. Considering the above bounds, we take / = 4 TeV or / = 5.6 
TeV with a large tan (3 to illustrate our results. 

In Fig. [31 we take several values of tan f3 and plot the hadronic cross section of Higgs-pair 
production at the LHC versus the Higgs boson mass. We find that compared with the SM 
prediction, the cross section in the SLHM can be significantly enhanced for a large tan f3. For 
example, with tan/?=18 (25) and / = 4 TeV (5.6 TeV), the cross section can be enhanced 
by 80% for = 110 GeV. Of course, for the perturbation to be valid, tan f3 cannot be too 
large for fixed /. As shown in Eq. f|T5|) . the correction to the Higgs VEV is proportional 
to tan 2 (3vq/ f 2 . If we require 0{v\j f^/O^v 2 / f 2 ) < 0.1 in the expansion of v, the value of 
tan/? should be below 20 (28) for / = 4 TeV (5.6 TeV). For a larger /, the value of tan/3 
can be larger and cancel partially the suppression of v/f. Therefore, the maximal value of 
the cross section does not always decrease with increasing of /. 
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IV. FINAL STATES OF HIGGS-PAIR PRODUCTION 



The Higgs-pair production can give various final states, depending on the decay modes 
of the Higgs boson. The SLHM corrections to the tree-level decays h — > //, WW, ZZ are 
mainly from the corresponding modified couplings: 

T{h - XX) = T{h - XX) SM (g hXX /g™ x ) 2 , (16) 

where XX denotes WW, ZZ or fermion pairs, and T(h — > XX) sm is the SM decay width. 
Qhxx and g^xx are the couplings of hXX in the SLHM and SM, respectively. The couplings 



ghww and ghzz can be found in 



12| 



For the low Higgs mass, the loop- induced decay h gg will be also important. In 
addition to the top quark loops, the loops of new heavy quarks (T, D, S) come into play. 
For another important loop-induced decay mode h —>■ 77, in addition to the contributions of 
top quark and W boson, the new charged heavy fermions (T, D, S) and gauge bosons W r± 
will make contributions. Following the approach in [lj]], the partial decay width of h — > 77 
can be calculated at one-loop level. For the SM decay channels, the relevant higher order 

n 

QCD and electroweak corrections are considered using the code Hdecay [20]. 

In addition to the SM decay modes, the Higgs boson in the SLHM has two new important 
decay modes, h — > rjr] and h —>■ Zrj, in the kinematically allowed parameter space. Their 
partial widths are given by 



r(h — > 7777) 



A' 2 v 2 








ml 




32vr/ 2 



2 / ml m 2 



nh^Z V ) = -^- A 3 / 2 1,3,-| , (17) 



ml m 2 h 



where x v = 4m 2 /ml and A(l, x, y) = (1 — x — y) 2 — 4xy. These two decay channels can be 



dominant in the allowed parameter space [12| and provide some new signatures of Higgs-pair 
production. 

Fig. H] shows the decay branching ratios of Higgs-pair versus the Higgs boson mass (we 
only plot the decay modes with branching ratio above 0.1). We see that the dominant decay 
channel is hh — > WWWW for 150 GeV < m h < 200 GeV, similar to the SM prediction; but 
for Higgs mass below 140 GeV, the decay hh — > r\r\r\r\ will dominate over hh — > bbbb which 
has the largest branching ratios in the SM. Besides, the decays hh — > 667777 and hh — > f]f]WW 
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FIG. 4: The decay branching ratios of Higgs-pair as a function of the Higgs boson mass. 
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FIG. 5: The rates of a(pp -> hh) x BR(hh -> bbr+T'), a(pp -> hh) x BR{hh -> 6677) and 
a(pp -> hh)xBR(hh -> WWWW) in LHT-I, LHT-II and SLHM, normalized to the SM prediction. 
Here, tan/3 is fixed as 18 and / is taken as its lower bound, which is 4 TeV for the SLHM and 500 



GeV for LHT-I and LHT-II 



2i. 



can also be important, whose branching ratios can be much larger than the decay hh — > bbbb 
in some part of the parameter space. 
The decays of rj have been studied in 



22 



23|. For 10 GeV < m v < 100 GeV, 77 decays 



mainly into bb, t + t~ or gg. The branching ratio of 77 — > t + t - is about 10% of 77 — > 66. With 
increasing of m^, the branching ratios of 77 —> 66 and 77 — > r + r~ decrease while the decay 
77 — > gg increases and may surpass the ratio of 77 — > 66. 

In the SM the promising channels are pp hh —>■ bbr + r~ (6677) for rrih < 140 GeV 



24j and pp -> /i/i -> WWWW for 150 GeV < m fc < 200 GeV [8|. In Fig. [5] we plot the 
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rates of a(pp — > /i/i) x BR(hh — > 66t + t~), cr(pp — > /i/i) x BRihh — > 6677) and crQop — > 
fo/t) x BR(hh — > WWWW) normalized to the SM predictions, and compare the SLHM 
results with the predictions of two types of littlest Higgs models with T-parity (LHT-I and 



LHT-II). The detailed descriptions of LHT-I and LHT-II can be found in [3J, |25J. We see 
that for rrih < 130 GeV all the three rates can be enhanced sizably in LHT-I and LHT-II, 
but suppressed significantly in the SLHM. For the larger value of rrih, both the SLHM and 
LHT-I/LHT-II can enhance sizably the SM predictions (in the SLHM for rrih > 150 GeV, 
while in the LHT for m h > 170 GeV). 



V. CONCLUSION 



In the framework of the simplest little Higgs model (SLHM), we studied the production of 
a pair of neutral CP-even Higgs bosons at the LHC and obtained the following observations: 
(i) The Higgs-pair production rate in the SLHM can be significantly larger than the SM 
prediction; (ii) For a low Higgs mass the dominant decay mode of Higgs-pair is hh — > 11111111 
(i] is a CP-odd scalar) while hh — > bbrji] and hh — > rjrjWW may also have sizable ratios; 
(iii) For a light Higgs boson all the rates of pp — >• hh — > bbr + T~, pp —>■ hh —>■ 6677 and 
pp — > hh — > WWWW can be sizably enhanced in the littlest Higgs models but severely 
suppressed in the SLHM; while for an intermediately heavy Higgs, all the three rates can be 
sizably enhanced in the littlest Higgs models and the SLHM. 
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